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ABSTRACT 


The  theory  of  Singular  Lines,  constant  damping  ratio-constant 
undamped  natural  frequency  lines,  is  derived. 

A  limitation  of  the  parameter  plane  method  for  characteristic 
polynomials  whose  coefficients  are  linear  functions  of  two  variable 
parameters  which  results  in  undetermined  roots  is  described.   The 
addition  of  singular  lines  to  the  parameter  plane  diagram  specifies 
these  roots  allowing  solution  for  all  roots  of  a  given  polynomial. 

A  general  method  of  solving  for  singular  lines  is  developed  and 
rules  for  predicting  the  existence  of  such  lines  are  stated.   A  computer 
program  which  solves  for  and  graphically  displays  singular  lines  in 
addition  to  constant  zeta,  omega,  and  sigma  loci  is  presented. 

Singular  lines  are  considered  in  terms  of  dominance  and  macro- 
scopic root  sensitivity.  A  dominant  root  line  in  the  parameter  plane 
is  illustrated. 

Examples  which  demonstrate  the  application  of  singular  line 
theory  to  linear,  multivariable  control  systems  are  presented. 
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I  INTRODUCTION 

Prior  to  1959  analysis  and  synthesis  of  control  systems  were 
carried  out  primarily  by  root  locus  and  frequency  response  techniques.  [1] 
Mitrovic's  Method  [2],  which  was  introduced  in  1959,  allowed  the  deter- 
mination of  a  system's  root  locations  as  two  coefficients-  of  the 
characteristic  equation  were  varied.   The  method  specified  the  variable 

coefficients,  B-  and  B  ,  as  the  two  lowest  order  terms  of  the  character- 
1      0 

istic  equation.   For  example: 

a  sn  +  a  ,sn-1  +  +  a„s  +  B.s  +  B  =0  (1-1) 

n      n-1  2      1      o 

where: 

a,      k  =  0,  1,  ...,n  are  real  constants 
k 

B-  and  B  are  real  variables 
1      0 

Application  of  the  method  resulted  in  constant  zeta,  omega,  and  sigma 

curves  plotted  on  the  B  -  B..  plane  specifying  the  roots  of  equation 

1-1  for  any  choice  of  C  and  oj   .   The  Coefficient  Plane  Method  includes 
J  n 

Mitrovic's  original  work  and  extension  of  his  method  to  include  variation 

of  any  two  of  the  coefficients  of  the  system  characteristic  equation. 

In  1964,  Siljak  [3]  extended  Mitrovic's  work  to  control  systems 

in  which  two  adjustable  parameters  appeared  linearly  in  the  coefficients 

of  the  characteristic  equation.   In  the  linear  system  case,  Siljak 

considered  the  characteristic  equation: 

n     , 
f(s)  =  I      a,  s*  =  0  (1-2) 

k=0  K 


where:     ak  =  bka  +  ck3  +  dk  (1-3) 

and:     b.  ,  c,  ,  d.     are  real  constants 
k   k   k 

a,  3  are  real,  variable  parameters 

This  method  produced  constant  zeta,  omega,  and  sigma  curves  plotted  on 
the  a  3  plane  specifying  the  roots  of  equation  (1-2)  for  any  choice  of 
a  and  3.   The  Parameter  Plane  Method,  as  Siljak  named  it,  is  a  simple 
procedure  for  factoring  characteristic  polynomials  and  displaying  the 
results  in  a  parameter  plane  diagram. 

The  Parameter  Plane  Method  was  extended  in  1965  to  a  general  case 
in  which  coefficients  of  the  characteristic  equation  are  a  nonlinear 
combination  of  two  adjustable  system  parameters.   [A],  [5],  [6]   Specif- 
ically, the  product  case  where  the  coefficients,  a,  ,  of  equation  (1-2) 


were: 


ak  -  bka  +  ck3  +  tytS  +  dk  (1-A) 


where: 


b.  ,c.  ,  h,  and  d.     are  real  constants 
k  k   k     k 

a,  3  are  real  variable  parameters 

A  further  extension  of  the  method  to  include  characteristic  equation 
coefficients  of  quadratic  form  is  limited  to  the  specific  case  of  third 
order  systems.   [7] 

The  Parameter  Plane  Method  was  conceived  as  an  approach  to  the 
analysis  and  synthesis  of  feedback  control  systems.   It  allows  the 
designer  to  obtain  information  about  system  relative  stability  and  the 
effect  of  parameter  adjustments  on  stability.  By  adjusting  pole-zero 
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locations  in  the  system  transfer  function  the  designer  maintains  control 
over  both  transient  and  frequency  responses.   [8]   The  method  is 
particularly  applicable  to  multiparameter,  multiloop  control  systems 
with  more  than  one  adjustable  parameter  since  it  is  based  on  the  study 
of  the  system  characteristic  equation  written  in  a  general  form.   In 
general  the  parameter  plane  method  may  be  applied  effectively  to  any 
engineering  problem  in  which  it  is  necessary  to  determine  how  variations 
of  parameters  in  the  characteristic  equation  effect  root  locations. 
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II  THE  PARAMETER  PLANE  METHOD 

This  section  will  review  Siljak's  derivation  of  the  parameter 
plane  method  for  the  linear  case  [3],  discuss  the  use  of  the  parameter 
plane  for  system  stability  analysis,  and  illustrate  the  primary  rules 
for  mapping  points  and  curves  from  the  s-plane  to  the  a  3  plane. 

Parameter  Plane  Equations 

Consider  the  characteristic  equation: 

f  (s)  =  I     a,  sK  =  0  (2-1) 

k=o   R 


where: 


the  coefficients  a   (k=o,l,  ...,  n)  are  real  and  s  is 
the  complex  frequency  variable: 


jw   ~Jl-;2  (2-2) 


n 


where  w  is  the  undamped  natural  frequency  and  C  is  the  damping 
n 


ratio 
By  forming  powers  of  s  it  can  be  shown  [3]  that: 


k 

s   =   0) 


where: 


nk  [Tk(-0  +  j  "^1-C2  Uk(-03  (2-3) 


Tk(-0  =  (-l)k  Tk(C)  (2-4) 

Uk(-5)  +  (-l)k+1Uk(0 

The  U,  and  T,  are  Chebyshev  functions  of  the  first  and  second 
kind  respectively.   They  are  given  by  the  recursion  formulae  [3]: 
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W°  - 2  V?)  +  Tk-i(?)  =  ° 


<w?)  - 2  v?)  +  \-l«>  -  ° 


(2-5) 


where: 


T0(C)  =  1  UQ(0  =  0 

T1(C)  =  5  ^(C)  =  -1 

Substituting  equation  (2-3)  in  equation  (2-2)  and  equating  the  real  and 
imaginary  parts  of  the  resulting  equation  to  zero  independently  gives: 

I     Vn  k  Tk(_°  =  ° 
k=0  R  n    k 


(2-6) 


I      Vn  k  Uk(_C)  =  ° 
k=0  R  n    k 


From  equations  (2-5)  we  obtain  the  interrelation  equation: 

TkU)  =  C  Uk(C)  "  Uk-l(C)  (2"7) 

Substitution  of  equations  (2-7)  and  (2-4)  in  equations  (2-6)  produces, 
after  simplification,  equations  in  one  Chebyshev  function: 


k        k 


X   ('X)  Vn   Uk-1(?>  "  ° 
k=0 


(2-8) 


I      (_1)k  Vn  k  V°  =  ° 
k=0       k  n    k 

Assume,  as  Siljak  did,  that  the  coefficients  a  of  the  character- 
istic equations  are  linear  combinations  of  parameters  a  and  $  of  the  form: 


ak  =  V  +  cke  +  dk  <2"9) 
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where 


b,  ,  c.  and  d.  are  real  constants 
k   k      k 


a  and  6  are  variable 


Substitution  of  equation  (2-9)  in  equation  (2-8)  yields: 


B  a  +  CJ  =  -D 


B2a  +  C23  =   -D2 


where: 


(2-10) 


Bl  =  X    (A"\lW  B2   =  I    (-1)kVn\^> 

k=0  k=0 


ci  -  X  H)V.Vi«  c2  •  Jn  <-»kv«V«      (2-u) 

k=0  k=0 


Dl  "  jn   (-1)kVn\-l(?)  D2  *  jn   t^\ft) 

k=0  k=0 


The  functional  dependence  of  B..  ,  B  ,  C  ,  C„,  D  ,  and  D_  on  z,   and  to  is 
omitted  in  equations  (2-10)  and  succeeding  equations  for  simplicity  of 
notation.   A  Table  of  the  Chebyshev  functions  of  the  second  kind,  U,  (O 
is  given  in  Appendix  I.   The  numerical  value  of  U,  (£)  for  selected  £  is 
given  in  Appendix  II. 

Application  of  Cramer's  Rule  for  solution  of  simultaneous  linear 
equations  to  equation  (2-10)  yields  the  desired  parameter  plane  solution 
equations : 

C1D2  -  D2D1  e  .  B2D1  -  B1D2  (2_12) 

B1C2  "  B2C1  B1C2  "  B2C1 
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As  defined  by  Siljak  [3]  the  parameter  plane  or  a$  plane  is  a 
rectangular  coordinate  plot  with  a  the  abscissa  and  $  the  ordinate. 
Equations  (2-12)  give  a  and  3  as  functions  of  C  and  to  .   Fixing  £  =  c; . 

and  varying  to   from  0  to  infinity  in  equations  (2-12)  produces  a  curve 

in  the  a3  plane,  a  constant  zeta  curve.   This  curve  specifies  the  a3 

pairs  which  will  cause  equation  (2-1)  to  have  a  pair  of  complex  roots 

with  the  required  damping  ratio  C, . .      Similarly,  fixing  to  .  and  varying 

£  from  -1  to  +1  in  equations  (2-12)  produces  a  curve  in  the  parameter 

plane,  a  constant  omega  curve,  specifying  the  a3  pairs  which  will  cause 

equations  (2-1)  to  have  a  pair  of  complex  roots  with  the  required 

natural  frequency  to  ..  Equations  (2-12)  thus  allow  mappings  of  points, 

excepting  real  axis  points,  from  the  s-plane  to  the  parameter  plane. 

For  real  axis  points  in  the  s-plane,  replace  s  in  equation  (2-1) 

by  s  =  -a  to  obtain: 

n       , 
f (s)  =  I     a  ,(-o)K  =  0  (2-13) 

k=0  K 

Now,  substitute  equation  (2-9)  for  a,  to  obtain: 

n  , 

I      (b,  a  +  c,  3  +  d,  )  (-a)K  =  0  (2-14) 

k=0   R     *     K 

Simplifying  equation  (2-14)  yields: 

aB  (a)  +  3C  (a)  +  D  (a)  =  0  (2-15) 

where : 


k 
a 


B(a)  =  I      (-l)k  b 
k=0 


n        b-  ^   L- 

C(a)  =  I      (-I)*  d,  a  (2-16) 

k=0       k 


k  .  k 
a 


D(a)  =  I      (-1)*  d  ■ 
k=0 
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For  a  given  value  of  a  the  functions  B(a),  C(a),  and  D(a)  are  constants 
and  equation  (2-15)  is  the  equation  of  a  straight  line  on  the  a3  plane, 
the  locus  of  a3  points  corresponding  to  real  roots  s  =  -a. 

Stability  Analysis  and  Mapping 

Absolute  stability  analysis  of  a  linear  control  system  consists 
of  determining  the  existence  of  roots  of  the  characteristic  equation  in 
the  right  half  of  the  s-plane.   Similarly,  relative  stability  analysis 
consists  of  determining  the  existence  of  roots  within  specified  areas  of 
the  s-plane,  generally  defined  by  constant  X,   or  constant  oj  loci. 
Mapping  of  the  specified'  areas  from  the  s-plane  to  the  aB  plane  using 
equations  (2-12)  and  (2-15)  allows  a  designer  to  choose  or  adjust  para- 
meter values  so  that  characteristic  equation  roots   lie  within  required 
areas  thereby  ensuring  stability  requirements  are  met.   The  rules  and 
graphical  techniques  for  the  mapping  process  and  the  properties  of  the 
mapping  are  discussed  in  detail  by  Siljak  [3]  and  Thaler  [9].   The  inter- 
pretations of  the  curve  on  the  parameter  plane  and  determination  of 
stability  is  not  a  simple  process;  the  reader  is  again  referred  to 
references  [3]  and  [9]  for  detailed  information. 

The  major  mapping  rules  will  be  illustrated  in  the  following 
example. 
Example  I : 

Consider  a  linear  system  whose  characteristic  equation  is  the 
fourth  order  polynomial: 

f(s)  =  0.04  s4  +  0.34  s3  +  (0.2a  +  1.12)s2  + 

(2-17) 
(0.5a  +  3  +  1.7)s  +  2g  +  1  -  0 
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From  the  parameter  plane  diagram  for  the  polynomial  shown  in  Figure  1 
we  can  make  the  following  observations: 

1)  An  s-plane  complex  point  maps  into  a  single  a3  plane  point. 

2)  An  s-plane  real  axis  point  maps  into  a  straight  line,  a 
constant  sigma  line,  on  the  a$  plane.   Real  roots  may  be 
evaluated  from  the  a  lines.   For  example  the  point  M(2. 0,0.83) 
on  the  a  =  1.5  line  has  a  root  s  =  -1.5. 

3)  A  constant  damping  ratio,  natural  frequency  or  settling  time 

contour  maps  into  a  constant  C»u>  or  Cw  curve  on  the  a3 
r  n      n 

plane.   Complex  roots  of  the  polynomial  for  specific  values 
of  a  and  3  are  determined  by  the  £  and  to  curves.   For 
example  point  M(2.0,  0.83)  is  at  the  intersection  of  the 
C  =  0.5  and  to  =3.0  curves,  therefore  the  polynomial  has 
complex  roots  at  s  =  -(0.5)  (3.0)  +  j  3.0~Wl-(0. 5)2 


=  1.5  +  j  2.6 

4)  Constant  sigma  lines  on  the  a3  plane  are  tangent  to  the 

C,   =  1  curve  at  w  =1.5. 
n 

5)  An  s-plane  stable  area  maps  into  an  area  of  the  a$  plane  in 
which  all  the  roots  of  the  polynomial  have  negative  real 
parts,  i.e.,  into  a  stable  root  area.   For  example,  consider 
the  area  in  the  left  half  s-plane  bounded  by  the  origin,  the 
radial  £  =  0.5  lines,  and  the  infinity  of  the  s-plane. 
Setting  s  =  0  in  equation  (2.17)  yields  6  =  -0.5  as  the  a   =  0 
boundary  in  the  a3  plane.   The  C  ■  0.5  radial  lines  map  into 
the  C  =  0.5  contour  of  Figure  1.   The  enclosed  area  on  the 
plane  of  Figure  1  is,  therefore,  the  desired  stable  root  area, 
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FIGURE  1. 

PARAMETER  PLANE  EXAMPLE 
4TH,  ORDER  POLYNOMIAL 


Ill  BASIS  OF  THIS  STUDY 

This  section  introduces  a  heretofore  unknown  limitation  of  the 
parameter  plane  method  which  results  in  undetermined  roots,  illustrates 
the  limitation  with  a  specific  example,  describes  initial  investigation, 
and  then  formulates  a  hypothesis  as  to  the  cause  of  undetermined  roots. 

Sections  IV  and  V  develop  the  mathematical  theory,  solution 
method,  and  computer  implementation  extending  the  parameter  plane  method 
to  ensure  determination  of  all  roots  of  a  given  polynomial.   Sections  VI 
and  VII  consider  this  extension  of  theory  in  terms  of  dominance  and  root 
sensitivity.   Section  VIII  contains  examples  of  analysis  of  specific 
linear  control  systems  which  indicate  the  applicability  and  potential 
of  the  theory  developed.   Section  IX  concludes  this  study  with  comments 
on  work  completed  and  recommendations  for  further  investigation  of  para- 
meter plane  theory  and  methods. 

Limitation  of  Parameter  Plane  Theory — Undetermined  Roots 

The  parameter  plane  method  derived  by  Siljak  [3]  and  the  computer 
program,  PARAM  A,  written  by  R.  M.  Nutting  [8]  purportedly  solved  the 
problem  of  determining  all  roots  of  a  given  polynomial,  characteristic 
equation,  in  terms  of  variable  parameters  a  and  g  displaying  the  results 
on  the  parameter  plane  as  constant  C ,  u)  ,  and  £00  curves  and  constant  a 
lines.   In  most  cases  the  existing  parameter  plane  method  does,  in  fact, 
solve  for  all  roots  of  a  given  polynomial,  and  in  all  cases  such  infor- 
mation as  it  does  produce  is  correct.   There  are,  however,  situations 
in  which  the  present  parameter  plane  method  does  not  provide  a  solution 
for  all  roots  of  a  polynomial  nor  correctly  predict   the  existence  of 
complex  roots  in  specific  areas  of  the  aB  plane.   Example  II  is  an 
illustration  of  such  a  situation. 
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Example  II 

Consider  a  linear  control  system  whose  characteristic  equation  is 

s6  +  80s5  +  (20a  +  1600)s4  +  840as3  + 

2  (3-1) 

(1600a  +  4003)s  +  16003s  +  16006  =  0 

Solution  by  the  parameter  plane  method  results  in  the  parameter  plane 
diagram  shown  in  Figures  2  and  3.   The  use  of  two  figures  was  necessi- 
tated by  the  complexity  of  the  constant  z,   curves.   Although  it  is  not 
obvious  in  Figures  2  and  3  the  constant  £  curves  all  go  to  the  infinity 
of  the  parameter  plane  as  to  increases.   Although  the  C  =  1  curve  has 
discontinuities  it  effectively  circles  the  origin  clockwise  and  then  goes 
to  infinity  in  the  second  quadrant  of  the  af3  plane.   The  computed  values 
of  a  and  3  for  the  constant  £  curves  and  the  use  of  smaller  graph  scale 
for  3  confirms  these  statements. 

The  area  of  absolute  stability,  all  roots  in  the  left  half  s-plane, 
is  bounded  in  the  a3  plane  by  the  L,   =  0  curve,  the  positive  a  axis  and 
infinity. 

In  this  example  the  present  parameter  plane  theory  correctly 
solves  for  all  roots  of  the  polynomial  in  the  complex  root  area  of  the 
a3  plane  in  Figure  2.   In  the  real  root  area  bounded  by  the  c,   =  1  curve 
and  the  stable  area  limits  the  present  parameter  plane  method  does  not 
determine  all  roots  of  the  polynomial;  only  two  or  four  real  roots  are 
determined  by  constant  a   line  intersections  and  no  information  is 
presented  for  the  remaining  roots.   For  example,  for  point  M   (40,  300) 
of  Figure  2  the  roots  are 
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*1  -  -   2-U  s3  =  ? 

s0  =  -37.88  s.  =  ? 

2  4 

s   =  ? 

Similarly,  for  point  M„  (20,  26.6)  the  roots  are 

*x   =  -  2.77  s  =  ? 

s  =  -  9.77  s,  =  ? 

2  6 

s3  =  -27.5 

s.  =  -39.3 

4 

Existence  of  Undetermined  Roots  in  Parameter  Plane  Solutions 

Further  investigation  of  the  situation  outlined  in  Example  II 
yielded  the  following  results.   In  that  part  of  the  parameter  plane  for 
which  present  theory  predicted  only  real  roots,  specific  a  -  8  points 
were  found  to  have  either  one  or  two  pairs  of  complex  roots.   In 
addition,  certain  points  had  the  same  complex  pair  in  common.   For 
example,  for  point  M   (40,  300)  of  Figure  2  the  complex  roots  are: 

s3/4  =  -  1.06  -  j  1.83 

s  .  =  -18.94  -  j  31.24 

For  points  M   (20,  26.6)  and  M   (80,   113)  the  complex  roots  are 

sc/,  =  -  0.35  -  j  1.16 
5/6  J 

The  results  were  obtained  by  substituting  specific  a(3  pairs  in  equation 

(3-1)  and  solving  for  the  roots  of  the  resulting  equation.   Moreover, 

joining  points  M  and  M_  by  a  straight  line  and  investigating  specific 

points  along  the  line  showed  that  all  such  points  had  a  common  pair  of 

complex  roots.   Thus  a  line  of  constant  C-constant  a>  was  located  on  the 
r  n 

parameter  plane. 
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As  previously  stated,  the  existing  parameter  plane  method  correctly 
solves  for  all  roots  of  a  given  polynomial  in  the  majority  of  cases  and  in 
all  cases  the  root  values  which  it  does  determine  are  correct.   It  is 
therefore  obvious  that  existing  theory  as  derived  [3]  is  correct  insofar 
as  it  goes,  but  it  overlooks  or  dismisses  some  special  case.   Reviewing 
the  derivation,  Cf.  ante  pp.  12-16,  it  was  apparent  that  the  use  of 
Cramer's  rule  in  solving  equations  (2-10)  presumed  that  the  coefficient 
matrix  was  non-singular,  i.e.,  that  the  value  of  the  determinant  which 
formed  the  denominator  in  relations  (2-12)  was  non  zero,  for  all  cases 
of  interest.   This  then  seemed  to  be  the  part  of  existing  theory  which 
could  lead  to  undetermined  roots,  i.e.,  incomplete  solution  of  the 
problem,  and  this  hypothesis  formed  the  basis  for  the  remainder  of  this 
study. 
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FIGURE  1 

PARAMETER  PLANE  EXAMPLE 
6TH.  ORDER  POLYNOMIAL 
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FIGURE  3 

PARAMETER  PLANE  EXAMPLE 
6TH.  ORDER  POLYNOMIAL 
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IV 


SINGULAR  SOLUTIONS  OF  LINEAR 
PARAMETER  PLANE  PROBLEMS 


In  Section  III  a  heretofore  unknown  limitation  of  the  parameter 

plane  method  resulting  in  undetermined  roots  was  described.   Example  II 

presented  a  specific  illustration  of  the  problem  of  undetermined  roots. 

It  was  shown  that  at  least  one  constant  £  -  constant  to  line  existed  and 

n 

that  it  provided  solutions  for  some  of  the  roots  which  existing  para- 
meter plane  theory  could  not  determine.   It  was  then  hypothesized  that 
the  assumption  that  use  of  Cramer's  Rule  for  solving  the  simultaneous 
linear  parameter  plane  equations,  equations  (2-10),  produced  the  complete 
solution  for  all  cases  of  interest  excepting  that  of  real  roots  was  in- 
correct.  This  section  extends  parameter  plane  theory  to  include  the 
complete  solution  of  equations  (2-10)  by  solving  for  presently  undeter- 
mined roots.   The  theory  of  constant  C  -  constant  and  to  lines  defined 

J  n 

SINGULAR  LINES  is  derived,  a  tractable  singular  line  solution  method  is 
developed,  and  general  rules  for  determining  the  existence  of  such  lines 
are  stated. 

I.   SINGULAR  LINE  THEORY 


Mathematical  Basis 

The  solution  of  a  system  of  linear  non-homogeneous  equations  by 
matrix  methods  is  well  known  and  is  contained  in  any  standard  textbook 
of  linear  algebra  [10]  or  matrix  textbook  [11].   Consider  the  system  of 
equations: 


all  al2   •'•   aln 


a21  a22 


l2n 


ml   m2        mn 


—        ■■ 

xl 

hi 

X2 

= 

h2 

• 

• 

X 

n 

h 
n 

_        _ 

(4-1) 


25 


or,  in  compact  notation, 

AX  =  H  (4-2) 

By  way  of  review  the  following  definitions  and  statements  are  given 
without  amplification  or  proof: 

1.  A  =  [a..]  is  the  coefficient  matrix. 

2.  [AH]  is  the  augmented  matrix. 

3.  A  system  which  has  a  solution  is  said  to  be  consistent.   A  consistent 
system  has  just  one  solution,  a  unique  solution,  or  infinitely  many 
solutions. 

4.  A  system  which  has  no  solution  is  said  to  be  inconsistent.   The 
problem  leading  to  such  a  system  of  equations  is  indeterminate. 

5.  In  the  consistent  case,  the  coefficient  matrix  A  and  the  augmented 
matrix  [AH]  have  the  same  rank.   In  the  inconsistent  case,  they  have 
different  ranks. 

6.  When  one  or  more  equations  of  a  system  can  be  derived  from  another 
by  multiplication  of  all  of  their  terms  by  a  constant,  the  equations 
are  dependent  or  equivalent. 

7.  The  matrix  A  is  called  non-singular  if  its  rank  r  =  n,  that  is,  if 
|a|  #  0.   Otherwise,  A  is  called  singular. 

8.  The  system  of  equations  has  a  unique  solution  provided  the  common 
rank  of  matrices  A  and  [AH]  is  equal  to  n,  the  number  of  variables, 
that  is  provided  |a|  ^  0.   Thus  the  system  has  a  unique  solution  if 
A  and  [AH]  have  the  same  rank  and  if  A  is  non-singular. 

9.  In  a  consistent  system  of  rank  r  <  n,  a  solution  can  be  obtained  for 
r  variables  in  terms  of  the  remaining  n  -  r  variables. 
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Derivation  of  the  Singular  Line  Solution  Method 

Consider  the  linear  parameter  plane  equations,  equations  (2-10), 
which  comprise  a  system  of  linear  non-homogeneous  equations.   Transposing 
terms  and  rewriting  the  equations  gives 


B  a  +  C±&  =  -   D1 
B2a  +  C23  -  -  D2 


(4-3) 


The  coefficient  and  augmented  matrices  for  equations  (4-3)  are 


A  = 


C;  a 


(4-4) 


AH  = 


Bi   ci   "D: 
b2   c2 


(4-5) 


The  present  parameter  plane  method  solves  equations  (4-3)  for  those 
cases  in  which  the  common  rank  of  A  and  [AH]  is  two  as  follows: 


a  = 


-Di    Di 

-D    C2 


Bl   Cl 
B2   C2 


Bl  "Dl 
B2  -D2 


Bl   Cl 
B2   C2 


C1D2  "  C2D1 


B1C2  "  B2C1 


B2D1  "  B1D2 


B1C2  "  B2C1 


(4-6) 
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There  remain,  however,  the  cases  in  which  equations  (4-3)  are 
linearly  dependent.   In  these  cases  matrix  A  is  singular  rank  one; 
therefore,  the  above  method  of  solution  is  not  applicable  since  the 
determinant  of  A,  |a| ,  equals  zero.   To  obtain  a  solution  we  must  first 
choose  a  C  -  (i)  pair  such  that  matrices  A  and  [AH]  have  common  rank,  one 
and  then  solve  one  of  equations  (4-3)  for  one  unknown,  6,  in  terms  of 
the  other,  a. 

For  matrices  A  and  [AH]  to  have  common  rank  one  of  the  following 
conditions  must  apply: 


B_ 


B1C2  "  B2C1  =  ° 


(4-7) 


-D, 


-D, 


C1D2  ~  C2D1  =  ° 


(4-8) 


Bl  "Dl 
B2  -D2 


B2D1  -  Bl°2  "  ° 


(4-9) 


To  choose  a  z,   -  to  pair  which  satisfies  these  conditions,  expand  the 


determinant 


Bi    ci 

B    C2 


in  terms  of  z,   and  to  ,  substitute  a  specific 

n'  r 


value  of  £  =  £  ,  equate  to  zero,  and  solve  the  resulting  polynomial  for 

co  .   Real  roots  of  the  polynomial  and  the  specific  value  of  C  -  C  are 

n  v     *  r  s 

then  substituted  in  equations  (4-8)  and  (4-9) ;  those  values  of  to  which 

satisfy  the  equations,  to   ,  are  the  values  for  which  matrices  A  and  [AH] 

s 
are  both  rank  one. 

Substitution  of  a  z,     -to   pair  in  equations  (4-3)  and  solution 

s 
of  the  resulting  linearly  dependent  equations  for  parameter  6  in  terms 

of  parameter  a  yields 
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°1        "I  D2        B2 

—  -  —  a  =  --r--  —  a  (4-10) 


where: 


k=0  s  k=0  s 


Cl  -     I      ™\\  Vl'V  C2  -  Jn    ^^  Vn  Vs> 

k=0  s  k=0  s 


k=0  s  k=0  s 


(4-11) 


Equation  (4-10)  is  the  desired  singular  line  solution.   It  is  the  equa- 
tion of  a  straight  line  on  the  aS  plane  with  slope  -  B../C.  and  3  axis 
intercept  -  D  /C  .   All  points  on  the  line  are  maps  of  the  s-plane  point 

(C  i  w   )  thus  it  is  a  constant  C  -  constant  a)  line.   Substitution  of 
s   n  n 

s 

any  values  of  a  and  6  corresponding  to  the  infinity  of  a-3  pairs  defined 
by  equation  (4-10)  in  the  characteristic  equation  (2-1)  results  in 
s-plane  roots: 


+ 

-e». 


-  Id)     ^ 

n    J  n    \ 

s      s   l 


I-/ 


The  derivation  leading  to  equation  (4-10)  is  correct,  but  it  is 

by  no  means  a  practical  or  appealing  solution  method  since  the  procedure 

for  solving  for  £-oj  pairs  which  result  in  singular  lines  (c,    ,  oi   )  is 

n  o    n 

s 

both  tedious  and  time  consuming.   The  presence  of  Chebyshev  functions 
in  equations  (2-11)  produces  polynomials  of  order  n  and  n-1  in  w  where 
n  is  the  order  of  the  given  characteristic  equation  (2-1) .   Expansion 
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of  equations  (4-7),  (4-8)  and'  (4-9),  which  in  itself  is  a  tedious 
process,  leads  to  a  polynomial  of  order  2n-l  in  go  .   For  complex  systems 
with  higher  order  characteristic  equations  the  resulting  polynomial  may 
be  extremely  difficult  to  solve. 

In  order  to  derive  a  more  tractable  solution  method  consider  a 
characteristic  equation  of  the  form: 


f(s)  =  I     (b.o  +  c,  3  +  d.  )sk  =  0 


(4-12) 


k=0 


Setting  n=i  for  i=2,3;  computing  B- ,  B~,  C  ,  and  C,  in  equations  (2-11); 

expanding  equation  (4-7);  and  grouping  terms  yields: 

For  n  =  2: 

B1C2-B2C1  =    [b0VblC0]u)n  "    [b0C2'b2C0]2^n2  +    ^WlV  =  ° 


For  n  =  3: 


B1C2-B2C1  = 


[b0CrblC0]wn  -    [b0C2-b2C0]2;Wn     + 
[b0C3-b3C0](A;2-1)wn3  +   [blC2-b2Cl]a)n3 
[blC3"b3Cl]2;wn4  +   [b2C3_b3C2]a)n5 


-  0 


Details  of  the  expansion  for  n=2,  3,  4  and  the  extension  to  the  general 
case  are  contained  in  Appendix  III.   The  result  for  the  general  case  is 
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*wv 


*  k-i.*r^ 


(<tft-  fcj  <*  ♦ *  [%>%."  i;,<:-K'u$  *2"j 


fk  c-U  .lo^"' 


or,  in  compact  notation: 


(4-13) 


LC0  -  B0C1  -  I      ce.  . 
12    2  1    . L .    ii 


(4-14) 


where: 


CE  =  [ce„.]  is  an  (nxn)  upper  triangular 
matrix  whose  elements  are 


'bi-ici+j-rbi+j-icx-i]  (-d^Vk21  3"2 


i+j  <  n+2 
i+j  mi  n+2 


(4-15) 


Equation  (4-14)  results  in  a  polynomial  in  z,   and  to  of  order 

2(n-l)  in  co  .   Substitution  of  a  specific  £  =  C   and  solution  for  real 
n  s 

values  of  ui  gives  the  s-plane  coordinates  for  which  the  coefficient 
matrix  is  singular,  equation  (4-4).   The  real  values  of  w  and  the 
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specific  5  =  £  which  satisfy  equations  (4-8)  and  (4-9)  are  the 
s 

coordinates  (c,    ,  u  )  that  ensure  matrix  A  is  singular  and  matrices  A 
s   n 
s 

and  [AH]  have  rank  one.  c,     and  go  are  substituted  in  equation  (4-10) 

s 

which  is  now  a  linear  equation  in  two  unknowns  a  and  3.   The  resulting 

equation  defines  a  line  of  constant  z,  -  constant  u  ,  a  singular  line, 

on  the  a3  plane. 

The  above  solution  method  appears  cumbersome  when  described  but, 
in  fact,  it  proved  to  be  very  simple  to  apply.   To  illustrate  the  solu- 
tion method  Examples  III  and  IV  are  given  below.   In  common  with  the 
present  parameter  plane  method,  the  complexity  of  computation  increases 
with  the  order  of  the  given  polynomial  and  computer  solution  is  preferable 
for  higher  order-systems. 
Example  III 

Consider  the  characteristic  equation  given  in  Example  I: 

f(s)  =  0.04s4  +  0.34s3  +  (0.2a  +  1.12)s2  +  (0.5a  +  3  +  1.7)s  +  23  +  1  ■  0 

(4-16) 

For  ease  of  coefficient  component  identification,  display  the  terms  as 
follows : 


k 

0 

1 

2 

3 

4 

bk 

- 

0.5 

0.2 

- 

- 

Ck 

2 

1 

- 

- 

- 

dk 

1 

1.7 

1.12 

0.34 

0.04 

Apply  equation  (4-14)  to  compute: 


B1C2~B2C1  "  "%  "  0'4(-25)wn2  -  °'2\ 
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Simplifying  and  equating  to  zero: 


0)  "  -  4cu  +5  =  0 
n       n 


Choosing  £     =0.5,  for  example,  and  solving  for  the  values  of  u)  for 
s  n 


which  matrix  A  is  singular: 


w   -  2£u)  +5  =  0 
n       n 


2  -  /4-20   ,  +  ._ 
«  =  ^ =  1  -  j2 


Since  there  are  no  real  values  of  u   the  matrix  is  non-singular  for  all 

n  ° 

values  of  a)  ,  no  constant  ^-constant  co  line  exists,  and  the  present 
n  n 

parameter  plane  method  correctly  solved  for  all  roots  of  the  given  poly- 
nomial. 
Example  IV 

Consider  the  characteristic  equation  given  in  Example  II: 
f(s)  =  s6  +  80s5  +  (20a  +  1600)s4  +  840as3  +  (1600a  +  4003)s2  + 


16003  +  1600  =  0 


(4-17) 


Display  the  terms  as : 


k 

0 

1 

2 

3 

4 

5 

bk 

- 

- 

1600 

840 

20 

- 

Ck 

1600 

1600 

400 

- 

- 

- 

dk 

- 

- 

- 

- 

1600 

80 
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Apply   equation    (4-14)    to   compute: 


B,Ca-B,C( 


-03«V  t|WJ«Hf  -  [iwo(Mo)](45x-l)^     +[lo(.Mj(8SVS)«4 


-['*oo2W       *  [iwo(M©)]2Sa£     -[i«o©Gw)](«*-0q,f 


-  [400(^OJ|  t£    +  [400 M]  *  I  «Jf 


Simplifying  and  equating  to  zero: 


^wn4  -  (8?2+19)oj  3  +  (16£3+160£)w  2  -  (336£2+76)u>  +  320?  =  0 

Choosing  £  =0.5,  for  example,  and  solving  by  synthetic  division  for  u 
gives: 


a)  4  -  42w  3  +  164u>  2  -  320to  +  320  =  0 
n       n       n       n 


a)  =  2.1115,  37.889,  and  a  complex  pair 


Thus  the  matrix  A  is  singular  for  £  =  0.5,  u  ■  2.1115  and  for  £  =.0.5, 

co  =  37.889. 
n 

Substituting  the  first  pair  in  equations  (4-8)  and  (4-9)  leads  to  the 

results  0  =  0  confirming  that  (?  ,o>  )  =  0.5,  2.1115  is  a  singular  point. 

s   n 
From  equations  (4-11) :  s 

Bj^  =  (-l)2(1600)(2.1115)2(l)  +  (-l)3(840)(2.1115)3[2(0.5)]  + 
(-1)4(20)(2.1115)4[4(0.25)-1] 


=  -774.06 
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C1  =    (-l)°(1600)(2.1115)°(-l)+(-l) (1600) (2.1115)  (0)  + 
(-1)2(400)(2.1115)2(1) 


=  M83.4 


T>1  =    (-l)4(1600)(2.1115)4[4(0.25)-l]+(-l)5(80)(2.1115)5[8(0.125)-2]+ 
(-1)6(2.1115)6[16(.5)4-12(.5)3+1] 


=   3268.8 


Substituting   these  values   in  equation    (4-10)   yields 


3268.8  774.06 

183.4  183.4     a 


17.81     +       4.223  a 

which  is  a  line  of  constant  ^-constant  to  (.5,  2.1115)  whose  3  axis 
intercept  is  17.81  and  whose  slope  is  4.22.   Similarly,  the  values 
C  ■  0.5,  03  =  37.889  satisfy  equations  (4-8)  and  (4-9)  and  define  a 
singular  line: 

8  =  5740  +  75.777a 

Mitrovic's  Method  and  the  Coefficient  Plane  Method 

Singular  line  theory  must  be  compatible  not  only  with  the  para- 
meter plane  method  but  also  with  its  special  cases,  Mitrovic's  method 
and  its  extension  to  the  coefficient  plane. 

The  coefficient  plane  method  applies  when  any  two  coefficients  of 
the  characteristic  polynomial,  equation  (2-1),  are  variables  A  and  A 
where  n  _>  p  >  q  _>  0.   The  solution  parallels  that  for  the  parameter 
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plane,  outlined  in  section  II,  up  to  and  including  equations  (2-8). 


Substitution  of  variable  coefficients  A  and  A  in  equations  (2-8)  yields 

P     q  y 


(-1)PA  u)PU    (0  +  (-DqA  o>qU   (?)  =  -  I  (-l)kaku)kUk_1(0 


n 

I 

k=0 

i  p.q 


(-l)PApooPUp(;)  +  (-l)qAA(0  =  I    (-l)kak0)kUk(O 

j4  p.q 


(4-18) 


Considering  these  as  a  system  of  linear,  non-homogeneous  equa- 
tions in  two  unknowns,  A  and  A  ,  the  coefficient  matrix  is 

P     q 


(-1)P03PU    (O      (-l)qo)qU    (O 
p-1  q-1 


(-l)P0)PUp(O 


which  simplified  is: 


(-i)Vu  (O 


(4-19) 


(-0)) 


(-W)P_qU 


p-1    U  _ 


q-1 


(-0))P_qU 


(4-20) 


Considering  that  p  >  q,  inspection  of  a  representative  table  of  Chebyshev 
functions  of  the  second  kind,  Appendix  II,  shows  that  the  determinant  of 
the  coefficient  matrix,  equation  (4-9),  cannot  equal  zero.   Thus,  since 
the  matrix  is  non-singular  except  for  the  trivial  case  u  =0,  equations 
(4-18)  cannot  be  linearly  dependent  and  the  solution  for  A  and  A  ,  if 
it  exists,  must  be  unique. 

In  order  to  test  the  applicability  of  singular  line  theory  to 
coefficient  plane  problems,  recast  them  as  parameter  plane  problems. 
All  possible  cases  are  covered  if  we  consider: 


36 


A   (or  A  )   =  f(a) 

p    q 

(4-21) 

A   (or  A  )   =   f(g)  or  f(a,g) 
q     P 

In  any  case,  since  only  two  coefficients  of  the  characteristic  equation, 

(2-1) ,  contain  the  variable  parameters  the  expansion  of  the  coefficient 

matrix,  equation  (4-13),  contains  one  and  only  one  term.   The  result  is 

that  singular  line  theory  always  produces  the  solution  for  the  trivial 

case,  a)   =0,  and  no  singular  lines  are  defined.   Singular  line  theory 

s 
is,  therefore,  compatible  with  the  coefficient  plane  and  Mitrovic's 

Method. 

II.  CONDITIONS  FOR  EXISTENCE  OF  SINGULAR  LINES 

In  order  for  singular  lines  to  exist  the  linear,  non-homogeneous, 
parameter  plane  solution  equations,  (4-3),  must  be  dependent „   The 
necessary  conditions  for  this  are  that  their  coefficient  matrix  must  be 
singular  and  their  coefficient  and  augmented  matrices  must  have  common 
rank  one. 

Ideally  an  existence  theorem  which  clearly  defines  the  conditions 
on  the  characteristic  polynomial  necessary  for  singular  lines  should  be 
stated.   This  cannot,  however,  be  done  since  the  coefficients  and 
constants  in  the  parameter  plane  solution  equations,  equation  (4-3),  are 
each  involved  functions  of  n+2  quantities  where  n  is  the  order  of  the 
characteristic  polynomial,  that  is,  functions  of  b  or  c,  or  d  and  z, 
and  03  .   Certain  general  rules  which  indicate  the  likelihood  of  the 
existence  of  singular  lines  can,  however,  be  formulated.   They  are 
1.   The  variable  parameters  a  and  3  must  be  contained  in  more  than  two 
of  the  coefficients  of  the  system  characteristic  equation. 
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2.  Equation  (4-13)  indicates  that  for  a  characteristic  equation  of 
order  n: 

a.  For  each  value  of  £  there  is  a  maximum  of  2(n-l)  values  of  u 

n 

which  produce  singularities,  excluding  the  trivial  case  of  to  =0. 

b.  For  each  value  of  to  there  is  a  maximum  of  (n-1)  values  of  £ 

n 

which  produce  singularities. 
There  are  an  infinity  of  possible  values  of  t,   for  0  <  C  Jj,  1»  each 
associated  with  a  maximum  of  2 (n-1)  values  of  to  which  produce 
singularities.   Thus,  if  the  conditions  for  singular  line  existence 
are  met,  there  is  a  continuum  of  such  lines  spanning  the  real  root 
section  of  the  parameter  plane  stable  area. 

3.  Coefficients  of  a  and  $  in  the  higher  powers  of  s  in  the  characteristic 
equation  have  the  most  effect  on  the  possible  existence  of  singular 
points. 

A.   If  in  the  characteristic  polynomial  b  =c  =0  for  all  k  >_  M,  the 

maximum  possible  number  of  values  of  to  per  constant  C  curve  which 
r  n  r 

produce  singular  points  is  M-l  excluding  the  trivial  case  of  to  =0. 
5.   If  in  the  parameter  plane  diagram  the  stable  area  is  clearly  split 
into  a  complex  root  area  in  which  all  the  constant  C  contours 
eventually  go  to  infinity  and  a  real  root  area  which  contains  only 
constant  o  lines,  the  possibility  of  existence  of  singular  lines 
is  high. 
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V    COMPUTER  PROGRAMMING  OF  THE  SINGULAR  LINE  SOLUTION  METHOD 

Appendix  IV  is  a  listing  of  the  digital  computer  program  developed 

in  the  course  of  this  study.   The  program,  PARAM  S,  is  comprised  of  the 

major  parts  of  PARAM  A,  a  program  for  solution  of  the  linear  parameter 

plane  problem  written  by  R.  M.  NUTTING  [4] ,  and  a  new  sub-program  to 

implement  the  singular  line  solution  method  developed  in  section  IV  of 

this  study.   Computation  and  graphing  of  constant  t,   curves,  constant  u> 

curves,  constant  a   lines  and  constant  £  -  constant  u)  lines,  singular 

lines,  is  provided.   The  many  options  of  quantities  computed,  curves 

plotted,  and  output  data  printed  and  the  format  of  necessary  input  data 

J   j 
are  explained  in  the  comment  sections  of  Appendix  IV, 

The  computation  procedure  used  in  the  singular  line  sub-program 

for  a  characteristic  polynomial  of  the  form: 

n     ,  " 
f(s)  =  I     a,  sK  =   0  (5-1) 

k=0  R 

where: 

ak  =  V  +  cke  +  dk 

and  for  a  specified  value  of  damping  ratio,  £,  is  as  follows: 

1.  Read  input  data  L,   and  the  coefficients  of  the  characteristic  poly- 
nomial b,  ,  c,  ,  and  d,  where  k  =  1 ,  2 ,  . .  .  ,  n+1 . 

2.  Form  the  matrix  CE  =  [ce„.l    .  ,  n 

ijJ    i=l, 2,.,., n 


j=l,2,.. . ,n 


whose  elements  are 


(b.c.  L.-b.  ,  .c.)  i+j  <  n+2 

0  i+j  .1  n+2 
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3.  Form  the  matrix  CF  -    [cf      ]        i«l,2,...,n 

j"l»2 n 

whose  elements  are 

(-1)±+1[V±(0]  i  -  j 

0  i  +   j 

CF  is  a  signed  diagonal  matrix  whose  elements  are  formed  from  the 
Chebyshev  functions  of  the  second  kind  for  the  chosen  value  of  £. 

4.  Generate  the  AA  matrix  where: 

AA  -  CE  x  CF 
Sum  the  elements  of  the  AA  matrix.   The  resulting  sum  is  the  expan- 
sion of  (B.Cj-B-C.,),  equation  (4-13).   The  summation  process  groups 

the  terms  in  powers  of  oj  and  stores  the  coefficients  of  u  in 
r         n  n 

ascending  order. 

5.  Solve  for  U)  using  the  quadratic  formula  or  subroutine  POLYRT  as 
appropriate.   POLYRT  is  a  FORTRAN  60  subroutine  capable  of  solving 
for  the  roots  of  a  polynomial,  order  n  where  2  <  n  <  100,  with 
arbitrary  complex  coefficients.   The  subroutine  uses  Lehmer's 
method  [12]  to  approximate  a  root  and  then  improves  on  it  with 
Newton's  method.   POLYRT  is  a  relatively  slow  but  accurate  subroutine 
which  requires  no  supervision  or  initial  guesses  of  root  values.   For 

a  polynomial  with  repeated  real  roots,  the  results  often  have 

-3 
erroneous  complex  parts  which  are  generally  less  than  2.0x10 

6.  Test  the  output  of  POLYRT  retaining  only  real  roots;  these,  together 
with  the  specified  value  of  ;,  define  the  singular  points.  If  there 
are  no  real  roots  print  "Matrix  Non  Singular",  and  halt. 
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7.  For  each  real  value  of  w  and  the  specified  z,   generate  B,  ,  B~ ,  C,  , 

n  12   1 

C  ,  D  ,  and  D  ,  equations  (2-11);  test  that  the  conditions  for 
dependent  solution,  equations  (4-7),  (4-8),  and  (4-9)  are  satisfied. 
If  B  =B  =C  =C_,  print  "System  of  Equations  indeterminate.   Rank 
Zero."  re  loop  for  next  value  of  w  or  half,  as  appropriate,   If  either 
(4-8)  or  (4-9)  is  not  satisfied,  print  "System  of  Equations  Incon- 
sistent. No  Solution."  re  loop  for  next  value  of  w  or  halt  as 

r  n 

appropriate. 

8.  If  conditions  are  satisfied  computer  3  for  values  of  a  determined  by 
X  scale  and  IYRIGHT  input  data,  equation  (4-10),  and  plot  a  graph  of 
the  singular  line.   As  determined  by  program  options  print: 

a.  £,  b)  ,  XAXIS  intercept,  slope,  third  parameter,  and 

b.  roots  of  the  characteristic  polynomial  at  selected  points  on 
the  singular  line. 

Subroutine  POLYRT  is  again  used  in  solving  for  the  roots.   The  printout 
of  roots  serves  both  as  a  check  on  the  correctness  of  the  plotted 
singular  line  and  as  an  indication  of  the  movement  of  other  roots  as 
the  singular  line  is  traversed. 
Example  V 

Consider  the  characteristic  polynomial: 

s4  +  (20a+42)s3  +  (20a+4003+161)s2  +  (2440a+1600e+3280)s  + 

(2360a+16003+1600)  =  0  (5-2) 

The  output  graph  from  program  PARAM  S,  the  parameter  plane  diagram,  for 

this  polynomial  is  shown  in  Figure  4.   It  is  comprised  of  constant  L, 

curves,  constant  u)  curves,  and  singular  lines,  constant  £-  constant  to 

lines,  for  selected  values  of  t,   and  u  . 

n 
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The  computed  output  data  associated  with  each  singular  line  Is 
contained  in  Appendix  V.  The  two  forms  of  computer  printout  of  singular 
line  information  provided  by  program  PARAM  S  are  illustrated,  Table  III 
contains  the  values  of  ;,  w  ,  XAXIS  intercept,  and  slope  for  each  singular 
line.  Table  IV  includes  the  roots  of  the  given  polynomial  at  selected 
points  on  a  representative  singular  line. 

The  constant  5  curves  are  computed  for  values  of  m  from  .02  to 
200.  Although  it  is  not  obvious  in  Figure  4,  computer  data  printouts 
and  smaller  scale  graphs  show  that  the  constant  %   contours  describe 
short  arcs  in  the  second  quadrant  of  the  ot@  plane,  jump  discontinuous!)? 
to  the  fourth  quadrant,  then  follow  smooth  curves  through  the  fourth 
and  third  quadrants  prior  to  proceeding  to  infinity  as  shown. 

The  singular  lines  in  the  stable  area,  effectively  bounded  by 

the  l,   -  0  contour  and  the  postiive  a  axis,  each  define  a  line  of 

constant  s  -  constant  w  and  are  coincident  with  constant  0  lines.  For 

n 

example,  singular  line  43  for  which  c   «0.5,  w  « 12. 723  is  a  line  of  0-0 

s 
pairs  for  which  the  given  polynomial  has  roots: 


S  -  -£  W    -  JW    Ml"?  2 

s  n    J  n   \   s 

a       a 


-  -6.3614  -  jll.00 
and  s  •  -0 

-  -0.8077 

As  shown  in  Figure  4  the  remaining  real  root  moves  from  §*7i.5  to  s—28.1 
as  a  is  increased  from  -5.0  to  0  and  the  root  continues  to  move  to  the 
left  in  the  s-plane  as  a   is  Increased  further.  The  remaining  singular 
lines  are  contained  in  the  unstable  region  of  the  «g  plane;  that  is,  one 
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or  more  roots  of  the  given  polynomial  has  a  positive  real  part  for  any 
a- 3  pair  in  this  region.  As  indicated  on  the  parameter  plane  diagram, 
these  singular  lines  are  grouped  in  three  distinct  sets.  Two  of  the 
three  sets  contain  singular  lines  in  the  a3  plane  for  which  a  specific 
pair  of  complex  roots,  the  singular  pair,  are  in  the  right  half  of  the 
s-plane.  The  other  set  contains  singular  lines  which  specify  the 
singular  complex  pair  and  one  real  root  in  the  left  half  of  the  s-plane. 
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FIGURE  4. 


SINGULAR  LINES  ON  THE  PARAMETER  PLANE 
OF  A  4TH.  ORDER  POLYNOMIAL 
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VI  DOMINANCE  OF  SINGULAR  ROOTS 

The  concept  of  dominant  mode  analysis  and  design  has  wide  applica- 
tion for  linear  systems.   Most  system  criteria  in  the  s-plane  such  as 
damping  ratio,  natural  undamped  frequency,  settling  time,  and  peak  over- 
shoot have  meaning  only  for  a  dominant  mode  system.   The  no  zero,  two 
complex  pole  system  model  is  accepted  as  the  standard  model  for  s-plane 
analysis  of  linear  systems.   The  model  characteristic  equation  is 

2  2 

s  +  2Cto  s  +  oo  (6-1) 

n     n 

for  which  the  poles  are  located  at : 

s  =  -Cw.  -  ju)n"yi-?  (6-2) 

When  possible  a  designer  attempts  to  build  a  system  such  that 
the  characteristic  equation  approximates  equation  (6-1)  by  choosing  one 
pair  of  complex  roots  to  satisfy  equation  (6-2)  and  forcing  the 
remaining  roots  to  be  far  to  the  left  in  the  s-plane  or  placing  zeros  to 
cancel  undesired  poles  *   The  chosen  roots  then  dominate  the  transient 
response.   Since  the  resulting  response  is  essentially  that  of  a  second 
order  system,  C  and  co  may  be  taken  directly  from  the  dominant  root 
location.   The  peak  overshoot,  settling  time,  bandwidth,  etc,  may  then 
be  read  directly  from  a  set  of  standard  second  order  system  graphs.   In 
addition,  a  third  order  system  model  with  no  zero,  two  complex  poles  and 
a  real  pole  in  the  left  half  s-plane  near  the  origin  may  be  considered 
to  approximate  the  standard  second  order  model  for  analysis  purposes* 

For  systems  with  higher  order  characteristic  equations  the 
approximate  analytical  method  described  in  reference  [13]  may  be  used  for 
dominant  mode  design.   Essentially,  the  approximate  method  consists  of 
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writing  the  system  characteristic  equation  in  terms  of  its  parameters 

and  factoring  the  characteristic  equation  to  derive  two  equations,  one 

second  order  in  s.   The  values  of  system  parameters  are  then  chosen  to 

ensure  that  the  roots  of  the  second  order  equation  are  dominant,  and 

the  roots  of  the  other  equation  are  far  to  the  left  in  the  complex 

plane. 

Having  designed  a  dominant  mode  system,  if  the  values  of  ;  and  w 

n 

pertaining  to  a  singular  line  on  the  parameter  plane  diagram  (c  ,  w  ) 

_______  s   I* 

+       I 2  s 

define  s-plane  roots  8  ■  -5  co   -  jw   "\/l""5  which  are  equal  to  or  near 

8        8   " 

the  chosen  dominant  roots,  then  there  is  a  point  and  possibly  many  points 
on  the  singular  line  which  define  values  of  parameters  a  and  $  which  en- 
sure dominance.   For  example,  singular  line  43  of  Figure  4,  Example  V, 

for  which  5  -  0.5  and  to   *  12.723,  is  a  dominant  root  line.   Table  IV 
s  n         ' 

8 

in  Appendix  V  shows  that  any  point  on  the  section  of  this  line  in  the 
stable  area  defines  values  of  parameters  a  and  $  for  which  the  character- 
istic equation  (5-2)  has  the  desired  dominant  roots,  a  real  root  near 
the  origin,  and  a  real  root  far  removed  in  the  left  half  of  the  s-plane. 

When  the  fortuitous  coincidence  of  chosen  dominant  roots  and 
singular  line  fixed  roots  occurs,  as  in  Example  V  for  singular  line  43, 
a  dominant  root  area  may  be  shown  in  the  a@  plane.   The  boundaries  of 
this  dominant  root  area,  which  are  determined  by  the  degree  of  dominance 
required,  are  divergent  singular  lines.   For  example  singular  lines  33 
and  53  of  Figure  4  bound  a  dominant  root  area  for  which  dominancy  is 
between  9*5  and  11.0. 
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VII  SENSITIVITY  ANALYSIS  OF  SINGULAR  LINES 

There  have  been  many  sensitivity  functions  defined  for  the  effects 
of  parameter  variations  on  the  transfer  functions  of  linear  control 
systems.   The  root  sensitivity  functions  developed  by  Kokotovic  and 
Siljak  [14]  provide  the  most  general  solution  to  the  sensitivity  problem 
for  small  parameter  variations  in  linear  control  systems,  and  are  most 
applicable  to  parameter  plane  methods.   In  addition  the  macroscopic  root 
sensitivity  functions  defined  by  F.  H.  Hollister  [5]  are  applicable  for 
large  parameter  variations. 

Sensitivity  Equations 

This  formulation  is  an  abbreviated  composite  version  of  the 
derivations  given  by  Kokotovic  and  Siljak  [14]  and  F.  H,  Hollister  [5], 
and  it  presents  only  the  salient  points.   Consider  the  characteristic 
polynomial: 

I     a  sk  =  0  (7-1) 

k=0  K 

where  the  n  coefficients  are  functions  of  the  system 

parameters  q   (r-1,2, . „ , ,m) . 
A  variation  of  system  parameters  q  will  change  the  coefficients  a 
resulting  in  a  change  in  the  n  root  locations.   In  order  to  evaluate 
the  change  in  the  root  locations,  the  real  and  complex  roots  are 
considered  separately. 

Consider  the  i   pair  of  complex  roots  of  equation  (7-1)  as: 

r.  ,._  =  -c.w   -  jco  -\[l-Z2  (7-2) 

i,x+l     l  n.   J  n.  \j         l 
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and  the  j   real  root  of  equation  (7-1)  as: 


'j  -  0]  (7-3) 

Macroscopic  sensitivities,  the  sensitivities  of  the  roots  (7-2)  and  (7-3) 
to  large  (finite)  changes  in  parameter  q  ,  are  defined  as: 

ci    A5i 
S 


i,r     Aqn     q^O.i^r 


co.    A  to 
Si,r     Aq^     q^O.tfr  (7-4) 


o     Ao 


j,r     Aqr     q1-0,i^r 

S.   ,  macroscopic  damping  sensitivity,  is  a  measure  of  the  change  in 
the  damping  ratio  r     of  the  complex  root  pair  r.  due  to  a  finite  change 

to , 

in  parameter  q  .  S    ,  macroscopic  natural  frequency  sensitivity,  is 
a  measure  of  the  change  in  the  undamped  natural  frequency  to   of  the 


complex  root  pair  r  due  to  a  finite  in  parameter  q  »  SA  ^     ,  macroscopic 

j 


real  root  sensitivity,  is  a  measure  of  the  change  in  the  real  root  r 


due  to  a  finite  change  in  parameter  q  . 

Microscopic  sensitivities,  the  sensitivities  of  the  roots  (7-2) 
and  (7^3)  to  small  (infinitesimal)  changes  in  parameter  q  ,  are  defined 


as: 


?i 


8i.r  "  Hl/3"r 


to 

ni 


•i.r  '  »»Bl»«r  (7-5) 
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The  evaluation  of  microscopic  sensitivities  equations, (7-5)  is 
described  in  detail  in  reference  [14]  and  is  reviewed  in  reference  [5] . 
In  essence  it  consists  of  writing  the  characteristic  polynomial  as  two 
equations  in  the  Chebyshev  functions  of  the  second  kind,  equations  (2-6); 
computing  the  partial  differentials  the  functions  with  respect   to  the 
parameter  q.,  £ . »  and  w.;  rearranging  the  results  to  get  two  simultaneous 
equations  in  3o)./3q  and  3w./8q  with  coefficients  and  constants  which 
are  summations  in  U,  ,  U,  ,  ,  U,   and  U     and  a,  ;  and,  solving  for  the 
desired  sensitivities. 
Macroscopic  Sensitivity  of  Singular  Lines 

The  sensitivity  of  singular  lines  to  finite  changes  in  parameters 
a  and  3  is  determined  by  applying  equations  (7-4)  using  values  of  a,  3, 
C,  and  u  obtained  from  the  parameter  plane  diagram. 
Example  VI 

Consider  the  characteristic  equation  of  Example  V  and  the  para- 
meter plane  diagram,  Figure  4.   Assume  an  operating  point  M   (25,  9.3) 
on  singular  line  23  for  which  C  =  0.3  and  u  =12.02.   To  move  to 
singular  line  33  involves  an  infinity  of  complex  root  sensitivities 
since  the  singular  complex  roots  are  represented  by  an  infinity  of  points 
on  line  33;  however,  the  following  are  considered  to  be  of  primary 
interest: 

1.   Movement  of  the  operating  point  on  a  line  normal  to  the  terminal 
singular  line.   In  this  case  movement  to  point  M~  (24.6,  12.7)  on 

line  33  for  which  C  ■  0.4  and  m     =12.37.   The  macroscopic  root 

n 

sensitivities  are 


49 


0.1 
24.6-25 


=  -0.25 


n 


0.35 


0.35 


24.6-25  "  -0.40 


=  -0.875 


0.1 


12.7-9.3 


=  +0.032 


0.35 


0 


12.7-9.3 


0.35 
3.4 


=  +0.103 


2.   Parameter  3  fixed.   In  this  case  the  operating  point  moves  to  point 
M„  (17.7,  9.3)  and  the  macroscopic  root  sensitivities  are 


si  =  s  * 


17°;19q  =  -0.013     S*n   =  %2f  =  _o.048 
17.7-25  a    -7.3 


3.   Parameter  a  fixed.   The  operating  point  moves  to  point  M,  (25,  12.9) 
and  the  macroscopic  root  sensitivities  are 


SC     =  S  n 


6  12.9-9.3 


=  +0.028  S„n  =  %~  =  0.097 


3        3.6 


The  sensitivities  computed  in  Example  VI  show  that  the  frequency 
of  the  singular  line  complex  roots  is  more  sensitive  than  the  damping 
ratio  to  finite  changes  in  both  parameters  singly  or  together.   Inspec- 
tion of  Figure  4  shows  that  this  result  is  true  for  all  singular  lines 
in  the  stable  area  since  they  are  uniformly  divergent  with  increasing 
positive  slope  for  larger  damping  ratio  and  smaller  natural  frequency. 
In  addition,  Figure  4  shows  that  the  sensitivity  of  singular  complex 
roots  to  parameter  variations  is  decreased  as  the  operating  values  of 
a,  3,  or  a  and  3  are  increased;  that  is,  as  the  singular  lines  diverge 
sensitivity  is  decreased. 
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Complex  root  sensitivities  of  the  singular  lines  could  be  used 
effectively  in  self  adaptive  control  systems.   For  example,  consider  a 
system  with  characteristic  equation  (5-2)  for  which  the  desired  operating 
point  is  M-(29.6,12.7)  on  singular  line  33  of  Figure  4.   Consider  that 
either  parameter  a  or  3  changes  from  its  design  value  due  to  aging, 
breakdown,  or  external  disturbance  resulting  in  changes  of  both  e,  and 
u  of  the  complex  roots  of  the  system.   As  shown  in  Example  VI,  natural 
frequency  is  more  sensitive  than  the  damping  ratio  to  finite  changes  in 
both  a  and  3.   A  sensor  at  the  output  could  measure  the  system  response 
to  a  standard  input,  compare  this  with  a  model  generated  output,  and 
use  the  difference  in  transient  frequencies  to  determine  the  change  in 
natural  damping  of  the  dominant  roots.   A  controller  could  then  determine 
and  apply  the  change  in  parameter  3  necessary  to  drive  the  operating 
point  back  to  singular  line  33.   Since  the  complex  roots  are  dominant, 
the  real  root  near  the  origin  is  fixed,  and  the  remaining  real  root  is 
far  removed  in  the  left  half  s-plane  for  any  point  on  the  singular  line, 
the  specific  point  of  return  to  line  33  is  not  significant.   Parameter 
3  is  therefore  chosen  as  the  correction  element  since 


t^Aa-O  >  (SaWo  and  (S3^Aa=0  >  WaH 
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VIII  APPLICATION  OF  SINGULAR  LINE  THEORY 

TO  LINEAR  CONTROL  SYSTEMS 


In  this  section  selected  linear  multivariable  control  systems  are 
considered  in  sufficient  detail  to  demonstrate  the  applicability  and 
potential  of  singular  line  theory.   The  analysis  of  each  system  is  not 
complete;  discussion  is  limited  primarily  to  the  singular  lines  on  the 
parameter  plane  diagram  and  their  interpretation. 

The  basic  system  to  be  analyzed  is  an  inertially  stabilized  space 
vehicle.   [13]   Specifically,  the  system  is  comprised  of  two,  cross 
coupled  channels  each  with  reference  and  disturbance  inputs,  plant,  lead 
compensator,  coupling  element,  and  unity  feedback  path.   Throughout  the 
succeeding  examples  complete  symmetry  of  plants,  compensators,  and 
coupling  coefficients  is  assumed.   That  is: 


Gx  -  G2  -  G 


G  =  G  =  G 


cl   c2   c 


all  a22 


a12=  a21 


Examination  of  the  system  flow  graph,  Figure  5,  yields  the 
characteristic  equation: 

A  =  l+2ainG  G+(a  2-a  2)(G  G)2  (8-2) 

11  C      11    Ll  C 

Example  VI 

Assume  plant  and  compensator  transfer  functions: 
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G  = 


K 


(8-3) 


Kl  (y?s+1} 


Substitution  of  equations  (8-3)  in  equation  (8-2)  and  definition  of  para- 

2    2      2 
meters  q1  =  2  3^1^  and  q„  =  (a..   -a  _  )  (KK  )   gives  the  characteristic 

equation  as: 


(yx)  s  +  2yxs   +  (1+yx  q,  )s  +  (x+yx)q  s  + 

2    2 
Cq-jT*  q2)s   +  2xq2s  +  q2  =  0 


(8-4) 


Using  the  approximate  analytical  methods  of  reference  [13] 
compensator  parameter  values  of  x  =  0.5  and  y  =  0.05  are  chosen  to 
provide  dominant  roots  with  t,   =   0.5  and  u  =1.   Consider  parameters  q. 
and  q„  as  variables  a  and  6  in  the  parameter  plane  coefficient  equation, 
equation  (2-9) .   The  characteristic  equation  then  becomes 


A  =  s6  +  80s5  +  (1600+20a)s4  +  840as3  + 


(1600a+4006)s  +  16006s  +  16006     =  0 


The  system  transfer  functions  are 


(8-5) 


C^s) 
R-^s) 


q2    q2 


4         /      q2   2\ 
-  42s    [  80+40-^  s  +  ] 


ql    ql 


O.lq. 


(8-6) 


R2(s) 


4    3  /      a12   \  2   a12      al2 
s  +42s  +  80-10-^q.,   3+40-^,3+40-^, 

I      all  X  /      all  X   ^l?1 


O.lq, 


(8-7) 


Cx(s) 
D1(s) 


^(s) 
D2(s) 


;4  +  80s3  +  I  1600+10q;L  s2  +  420q  s  +  800q 


s  +  425s  +  80 


al2K  Kl 

2A 


(8-8) 


(8-9) 
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The  numerators  of  equations  (8-6) ,  (8-7) ,  and  (8-9)  are  real  functions 
of  one  relatively  simple  parameter  formed  from  the  gains  and  coupling 
coefficients  of  the  system,  i.e.,   2/q  ,   124l/a11,  and  12   1/2.   The 
zeros  of  these  transfer  functions  can,  therefore,  be  obtained  by  root 
locus  methods.   The  numerator  of  equation  (8-6)  is  a  function  of  para- 
meters K  and  q. ;  the  zeros  of  this  transfer  function  can  be  specified 
by  two  methods : 

1.  Fix  K  and  vary  q..  to  determine  the  root  locus. 

2.  Designate  q  =  a,  —  =  $  and  use  the  quadratic  coefficients  parameter 

i       K 

plane  method  to  obtain  a  parameter  plane  diagram. 
These  transfer  functions  are  included  simply  to  indicate  a  possible 
analysis  method;  no  discussion  of  them  will  be  given  since  it  is  outside 
the  scope  of  this  study. 

Figure  6  is  the  complete  parameter  plane  diagram  of  equation  (8-5) 
for  selected  values  of  £  and  w  .   Associated  computer  output  data  is 
contained  in  Appendix  V,  Table  V.   The  constant  £  and  constant  co  curves 
in  Figure  6  were  predicted  by  previous  parameter  plane  theory,  but  the 
singular  lines  were  not. 

For  system  specifications  requiring  an  underdamped  or  critically 

damped  response,  previous  parameter  plane  theory  restricted  values  of 

parameters  a  and  3  to  the  section  of  the  parameter  plane  bounded  by  the 

a  =  0  contour  and  the  C  =  1  contour.   Thus  values  of  a  and  6  which 

specified  an  operating  point  in  the  stable  real  root  area  were  prohibited, 

This  example  shows,  however,  that  the  singular  lines  define  operating 

points  at  which  complex  roots  exist  throughout  the  area  previously 

considered  solely  a  real  root  area.   Each  singular  line  is  tangent  to  an 

associated  constant  t  curve.   The  value  of  oj  on  the  constant  L,   curve  at 

n 
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the  point  of  tangency  defines  the  constant  L,  and  constant oj  of  the 
singular  line.  In  addition,  in  this  particular  case,  the  parameter 
plane  diagram  shows  that : 

1.  Each  singular  line  is  not  only  a  line  of  constant  c,   -   constant  oj 

n 

for  one  pair  of  complex  roots,  but  is  also  a  line  for  which  the 
value  of  one  real  root  is  constant.   Thus  the  singular  lines  are 
coincident  with  constant  a   lines.   Specification  of  the  two  system 
parameters  q.  and  q»  such  that  they  define  any  point  on  a  singular 
line  fixes  three  of  the  six  roots  of  equation  (8-5) .   The  remaining 
three  roots  vary  with  the  position  along  the  singular  line.   For 
example,  singular  line  101  for  which  £  =  0.9  and  co  =3.83  has 
fixed  roots: 

s±  2  =   -3.448  -  j  1.670 

s3   =  -33.103 

2.  The  intersection  of  two  singular  lines  defines  the  values  of  the 
system  variable  parameters  qn  and  q„  necessary  to  fix  all  six  roots 
of  equation  (8-5).   For  example,  operating  point  M  (64.5,  650)  which 
is  at  the  intersection  of  singular  lines  72  and  101  results  in 
characteristic  equation  roots: 

Sl/2  =  ~18'92  "  J25-22 

s„   =  -21.648  determined  by  line  72  and, 

*tu   =  "3.448  -  jl.67 
5/6  J 

s,   =  -33.103  determined  by  line  101. 

6 

This  was  to  be  expected  since  such  an  intersection  defines  unique 
values  of  a  and  6  so  that  the  coefficients  of  equation  (8-5)  are 
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all  constants.   Previous  theory  predicted  the  intersection  of  six 
constant  a  lines  rather  than  two  singular  lines  each  coincident  with 
a  constant  a   line  and,  in  fact,  six  intersecting  constant  a  lines  do 
not  exist. 

3.  Singular  line  41  for  which  C,   -   0.3  and  u  =1.26  clearly  defines  a 

n 

dominant  root  line  in  the  stable  area  of  the  parameter  plane. 
Despite  variations  in  parameter  a  and  6  along  the  line,  the  roots 
s  =  -0.3793  -  j 1.2062  remain  dominant,  the  real  root  o  =  -39.24 
remains  fixed,  and  the  other  roots  move  to  satisfy  the  characteristic 
equation.   Inspection  of  the  computer  data  printout  for  singular  line 
41,  Table  VI  of  Appendix  V  shows  that  the  singular  roots  have  a 
dominancy  of  about  three  at  point  M~(10,13.2),  but  as  the  operating 
point  moves  along  the  line  to  M  (80,  123)  the  dominancy  is  reduced 
to  approximately  one  point  six. 

4.  Theoretically  a  dominant  root  area  bounded  by  chosen  singular  lines 
diverging  from  singular  line  31  could  be  defined.   The  choice  of 
specific  lines  would  be  determined  by  the  degree  of  dominance 
acceptable;  for  example,  the  singular  roots  have  dominancy  of  between 
2.3  and  1.2  as  singular  line  31  is  traversed.   This  example  shows 
that  the  approximate  analytical  method  of  reference  13  is  at  best 
very  approximate.   The  required  dominant  roots  were  t,  =  0.5  and  oj  =  1 
with  a  dominancy  of  about  ten,  but,  in  fact,  on  singular  line  51  with 
C  =  0.5  and  w  =  1.055  the  dominancy  is  between  1,5  and  1. 

Figure  7,  Sensitivity  of  Singular  Lines  to  Parameter  Variations, 
is  the  parameter  plane  diagram  of  equation  (8-5)  for  values  of  c,   between 
0.5  and  0.6  and  £  between  0.8  and  0.9.   Inspection  of  the  diagram  yields 
the  following  information.   Maximum  sensitivities  of  c,   and  w  of  the 
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SINGULAR  LIMES  ON  THE  PARAMETER  PLANE 
OF  A  MULTIUARIABLE  ,   COUPLED  SYSTEM 
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singular  complex  roots  to  finite  variations  in  a,  3,  and  a  and  3  occur 
at  the  points  of  tangency  of  the  singular  lines  with  associated  constant 
z,   contours.   That  is,  points  M   (8.6,22)  for  £  between  0.5  and  0.6  and 
M   (2.175,  119)  for  £  between  0.8  and  0.9.   As  singular  lines  diverge 
from  these  points  the  macroscopic  singular  complex  root  sensitivities, 
equations  (7-4) ,  decrease. 

Consider  operating  point  M  (36.35,  135.6)  on  singular  line  11 
for  which  z,   =  0.5  and  u  =  2.1115.   Moving  the  operating  point  to 
M,(36.35,  141.2)  on  singular  line  21,  ?  -  0.51,  u  =  2.154  yields 
sensitivities: 

/s\       0.51-0.50     0.01 

lS3       "  141.2-135.6  "  5.6    °-0018 

V    Aa=0 

Aa=0 

That  is,  the  natural  frequency  is  approximately  four  times  more  sensitive 
to  changes  in  3  with  a  constant  than  the  damping  ratio.   Moving  the 
operating  point  to  M  (35.3,135.6)  on  singular  line  21  yields  sensitivities 

la/        35.3-36.35   -1.05 

A  3=0 

L\\  =o^i|  ^ 

I  a  I        -1.05 
N   7A3=0 

That  is,  the  natural  frequency  is  approximately  four  times  more  sensitive 
than  the  damping  ratio  to  changes  in  a  with  3  constant.   Moreover,  a 
finite  change  in  a  has  approximately  five  times  the  effect  of  an  equal 
change  in  3  on  the  natural  frequency  and  damping  ratio  of  the  singular 
lines. 
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Example  VIII 

Consider  that  the  specification  for  the  basic  system's  response 
is  changed  from  underdamped  with  one  overshoot  to  critical  damping. 
Assume  the  same  system  structure  and  transfer  functions.   The  system 
characteristic  equation  is  equation  (8-4) .   The  approximate  analytical 
methods  of  reference  13  now  require  compensator  parameters  f  =  0.707  and 
Y  =  0.05  to  place  dominant  roots  at  C  =  0.707  and  co  =  1.   Considering 
parameters  q1  and  q„  defined  in  Example  VII  as  variables  a  and   $,  the 
system  characteristic  equation  is 


s6  +  56.56s5  +  (20a+800)s4  +  593. 9as3  + 


(800a+4003)s2  +  1131.56s  +  8003    =  0  (8-10) 


The  parameter  plane  diagram  for  equation  (8-10)  is  shown  in 
Figure  8.   Computed  data  for  the  singular  lines  is  contained  in  Appendix 
V,  Table  VII.   The  constant  t,   curves  are  similar  to  those  in  Figure  6 
or  Figures  2  and  3  in  that  they  divide  the  stable  region  of  the  parameter 
plane  into  two  distinct  areas.   The  stable  region  is  bounded  by  the  C  =  0 
contour  and  the  positive  a  axis.   The  section  of  the  stable  region 
containing  the  intersecting  constant  c;  curves  was  previously  considered 
the  complex  root  region;  existing  parameter  plane  methods  correctly 
solved  for  all  six  roots  of  equation  (8-10)  in  this  region.   The  section 
of  the  stable  region  bounded  by  the  £  =  1  contour  and  the  positive  a  axis 
was  considered  the  stable  real  root  area.   Singular  lines  exist  in  this 
area.   Previous  parameter  plane  methods  determined  either  four  or  two 
real  roots  for  each  point  in  the  stable  real  root  area,  displaying  the 
reaults  as  intersections  of  constant  a   lines.   Inclusion  of  singular 
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lines  completes  the  parameter  plane  diagram  allowing  graphical  solution 
for  all  roots  of  equation  (8-10). 

Figure  8,  the  parameter  plane  diagram,  shows  that: 

1.  The  value  of  constant  w  for  a  singular  line  is  determined  by  the 

value  of  w  on  the  constant  £  curve  at  the  point  of  tangency  of  the 

singular  line.  For  example,  singular  line  51  (£-0.5,  w  -1.4929)  is 

n 

tangent  to  the  c-0.5  curve  at  point  M.(3.95,  3.91)  whereas  line  52 
is  tangent  to  the  ;  -  0.5  curve  at  point  M2(76.4,  1460.), 

2.  Each  singular  line  is  coincident  with  a  constant  a   line.  For 
example,  line  71  for  which  5-0.7  and  w  -2.0972  has  fixed  roots: 

8;L  2  -  -1.468  -  j  1.4977 
s3   -  -25.34 

3.  As  in  Example  VII,  the  intersection  of  two  singular  lines  define  the 
value  of  parameters  a  and  g  necessary  to  fix  all  roots  of  equation 
(8-10). 

4.  The  singular  complex  roots  become  dominant  at  singular  line  61,  that 
is  for  £-0.6.  The  amount  of  dominancy  is  increased  as  the  operating 
point  moves  to  singular  lines  71,  81,  and  91,  i.e.,  as  5  is 
Increased  to  0.9,  but  the  maximum  achieved  is  approximately  1.9. 

5.  The  natural  frequencies  and  damping  ratios  of  singular  lines  11 
through  101  are  more  sensitive  to  finite  changes  in  parameter  S 

than  to  changes  in  a.  Since  these  singular  lines  diverge,  effectively 
from  the  origin  of  Figure  8,  the  singular  root  sensitivities  increase 
with  respect  to  a   but  decrease  with  respect  to  I  as  the  operating 
point  6  value  increases  with  a  held  constant. 
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Example  IX 

Consider  the  basic  coupled  system  structure  of  Figure  5  but 
assume  plant  transfer  functions: 

G  -  — y* (8-11) 

(s>2s+2) 

Assume  the  compensator  transfer  functions  and  parameter  values  of 

Example  VII.   In  addition,  consider  system  parameters  q.  and  q«  as 

defined  in  Example  VII  as  parameter  plane  variables  a  and  $ .   The 

resulting  characteristic  equation  is 

s6  +  84s5  +  (20a+1928)s4  +  (880a+7048)s3  +  (3320a+400g+13444)s2  + 

(488Oa+16OO0+1312O)s  +  (3200a+1600@+1600)   -   0         (8-12) 

Figure  9  is  the  parameter  plane  diagram  of  equation  (7-12)  for 
selected  values  of  £  and  u  ,     Table  VIII  of  Appendix  V,  contains  the 
computed  data  for  the  singular  lines.  As  w  varies  from  0  to  *  the 
constant  £  contours  curl  around  the  origin  through  the  second,  third  and 
fourth  quadrant  of  the  a@  plane,  reverse  their  direction  producing 
sharp  spikes  in  the  lower  section  of  the  first  quadrant,  and  then  proceed 
to  infinity  via  the  first  quadrant  of  the  ag  plane.  For  clarity  only  the 
final  sections  of  these  contours  have  been  shown  in  Figure  9;  the 
sections  not  shown  are  similar  to  the  spikes  of  the  %   -  1  contour.  The 
x,   -  1  contour  is  a  particularly  confusing  contour  even  with  small  scale 
graphs  and  complete  data  printouts  available.  For  values  of  w  from  0.2 
to  2.06  it  describes  short  arcs  in  all  four  quadrants  of  the  a$   plane 
with  Intervening  discontinuous  jumps.   For  values  of  m     greater  than 
2.06  the  curve  is  as  shown  in  Figure  9  crossing  the  origin  and  proceeding 
to  infinity  via  the  second  quadrant. 
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The  bounds  of  the  stable  root  area  are  the  x,   =  0  contour  and  a 
line  determined  by  setting  s=0  in  equation  (8-12),  that  is: 

3200a  +  3  +  1600  =  0 

3  =  -(2a  +  1)  (8-13) 

Within  this  stable  area  there  are  four  distinct  sections  corresponding 
to  different  real  and  complex  root  combinations.   Previous  parameter 
plane  methods  solved  for  three  pairs  of  complex  roots,  two  real  roots, 
four  real  roots,  or  two  real  roots  and  a  complex  pair  in  these  sections 
of  the  stable  area.   The  inclusion  of  singular  lines  in  the  parameter 
plane  diagram  defines  the  remaining  roots  and  allows  graphical  solution 
for  all  six  roots  of  equation  (8-12)  at  any  point  in  the  stable  area. 
The  parameter  plane  diagram  and  associated  singular  line  data 
show  that : 

1.  Singular  line  values  are  accurately  determined  for  specified  values 

of  C  and  w  .   For  example,  line  51  for  which  C  =  0.6,  w  =  33.272; 
n  ft  n 

the  complex  roots  s  =  -19.963  -  j  26.618  are  correct  to  three 
decimal  places. 

2.  No  singular  line  can  be  considered  as  a  dominant  root  line  and 
therefore  a  dominant  root  area  cannot  be  defined. 

3.  The  singular  lines  are  technically  not  coincident  with  constant  a 

lines.   One  spike  of  the  5  =  1  contour  for  co  from  1.93  to  2.26 
r  n 

produces  a  mesh  of  double-valued,  constant  a  lines  which  encompasses 
the  stable  real  root  area.   The  result  is  that,  effectively,  each 
singular  line  is  coincident  with  a  double-valued  constant  a   line. 

4.  In  part  of  the  stable  area  singular  line  intersections  occur  solving 
for  four  roots  which  were  previously  undetermined. 
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5.  When  a  singular  line  is  extended  into  the  fourth  quadrant  of  the  a8 
plane  below  the  stable  area  boundary  defined  by  equation  (8-13)  the 
complex  roots  defined  by  the  singular  line  remain  fixed  but  the 
other  roots  shift  to  produce  a  right  half  s-plane  root.   For 
example,  since  line  73  crosses  the  stable  area  boundary  as  a  is 
decreased  from  40.0  to  20.0  a  system  root  moves  from  a  -  -2.02  to 
a  ■  +4.97. 

Example  X 

Consider  the  system  structure,  parameter  values,  and  parameter 

plane  variables  of  Example  VI  but  now  assume  plant  transfer  functions: 

G  -  ufir  (8-14) 

The  resulting  characteristic  equation  is 

s4  +  (20a+42)s3  +  (20a+400g+161)s2  +  (244Oa+16OO0+328O)s  + 

(1600a+16008+1600)  -  0  (8-15) 

The  parameter  plane  diagram  for  equation  (8-14),  Figure  10,  shows 
that  no  singular  lines  exist.  Values  of  £  between  0  and  1  were  Investi- 
gated. For  all  values  of  c  less  than  1  the  coefficient  matrix  of  the 
parameter  plane  solution  equations,  equation  (4-4),  is  non-singular. 

For  £»1  there  are  two  values  of  w  for  which  the  coefficient  matrix  is 

n 

singular,  but  for  both  values  the  parameter  plane  solution  equations 
(4-3)  are  inconsistent. 
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IX  CONCLUSIONS  AND  RECOMMENDATIONS 

FOR  FURTHER  WORK 


Singular  line  theory  is  a  small  but  meaningful  addition  to  the 
parameter  plane  method.   Singular  lines,  constant  damping  ratio  - 
constant  undamped  natural  frequency  lines,  occur  in  only  a  limited 
number  of  cases,  but  where  applicable  their  addition  to  the  parameter 
plane  diagram  defines  previously  undetermined  roots.   As  a  result, 
complete  solution  of  all  characteristic  polynomials  in  which  two 
variable  parameters  appear  linearly  in  the  coefficients  is  now  possible. 

The  conditions  for  existence  of  singular  lines  are  complex. 
General  rules  for  predicting  singular  line  existence  were  developed, 
but  an  existence  theorem  which  specifies  the  exact  combination  of 
coefficients  that  allows  singular  lines  should,  if  possible,  be  developed, 
Such  an  existence  theorem  would  facilitate  synthesis  and  design  of 
singular  line  systems. 

The  computer  program  developed  during  the  course  of  this  study, 
although  relatively  slow,  provides  an  accurate  means  of  solving  for  and 
graphically  displaying  singular  lines  in  the  parameter  plane  diagram  in 
conjunction  with  constant  £,  u  ,  and  a  curves.   More  efficient  program- 
ming techniques  and  the  use  of  a  faster  polynomial  root  solving  sub- 
routine should  be  considered. 

Dominance  of  singular  roots  results  in  a  dominant  root  line  on 
the  parameter  plane  diagram.   This  allows  choice  of  values  of  parameters 
a  and  (3  from  the  infinity  of  points  specified  by  a  singular  line  while 
simultaneously  ensuring  dominant  mode  operation  of  the  system. 

Sensitivities  of  the  damping  ratio  and  undamped  natural  frequency 
of  singular  roots  to  finite  changes  in  parameters  a  and  3  were  computed 
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for  specific  linear,  multi variable  control  systems.   It  is  considered 
that  measuring  the  movement  of  singular  roots  resulting  from  finite 
parameter  variations  could  prove  to  be  a  useful  concept  in  self-adaptive 
control  systems.   Since  movement  along  a  dominant  root  singular  line 
implies  little  change  in  operating  conditions,  the  adaptive  controller 
need  only  drive  the  system  back  to  the  singular  line  not  to  a  specific 
operating  point. 

Singular  line  theory  is  based  on  the  special  case  of  linearly 
dependent  parameter  plane  solution  equations.   The  theory  may  not, 
therefore,  be  extended  to  include  polynomials  whose  coefficients  are 
linear  functions  of  two  parameters  and  their  product. 


70 


X.  BIBLIOGRAPHY 


1.  Thaler,  G.  J.  and  Brown,  R.  G.   Analysis  and  Design  of  Feedback 
Control  Systems.   McGraw-Hill  Book  Co.,  Inc.,  New  York,  1960. 

2.  Mitrovic,  D.   Graphical  Analysis  and  Synthesis  of  Feedback  Control 
Systems,  I  -  Theory,  II  -  Synthesis.   AIEEE  Transactions,  Pt.  II 
(Application  and  Industry),  Vol.  77,  1958,  pp.  476-496. 

3.  Siljak,  D.  D.   Anaylsis  and  Synthesis  of  Feedback  Control  Systems 
in  the  Parameter  Plane,  I  -  Linear  Continuous  Systems.   IEEE 
Transactions  on  Applications  and  Industry,  Vol.  83,  No.  75,  November, 
1964,  pp.  449-458. 

4.  Siljak,  D.  D.   Generalization  of  the  Parameter  Plane  Method.   IEEE 
Transactions  on  Automatic  Control,  Vol.  11,  No.  1,  January,  1966, 
pp.  65-70. 

5.  Hollister,  F.  H.  Network  Analysis  and  Design  By  Parameter  Plane 
Techniques .   Ph.D.  in  EE  Thesis,  Naval  Postgraduate  School,  1965. 

6.  Miller,  A.  R.   Parameter  Plane  Analysis  of  Sampled  Data  Systems 
(With  Extensions  for  Continuous  System  Applicability).   Ph.D.  in  EE 
Thesis,  Naval  Postgraduate  School,  1965. 

7.  Paine,  J.  A.   Application  of  Parameter  Plane  Techniques  to  Control 
System  Compensation  (Identical  Double  Section  Cascade  Compensation) . 
M.S.  in  EE  Thesis,  Naval  Postgraduate  School,  May  1966. 

8.  Nutting,  R.  M.   Parameter  Plane  Techniques  for  Feedback  Control 
Systems.   M.S.  in  EE  Thesis,  Naval  Postgraduate  School,  1965. 

9.  Thaler,  G.  J.   Algebraic  Methods  for  Dynamic  Systems.   Unpublished 
notes  for  use  at  Naval  Postgraduate  School,  1965. 

10.  Murdoch,  D.  C.   Linear  Algebra  for  Undergraduates.   John  Wiley  and 
Sons,  Inc.,  New  York,  1957. 

11.  Ay res,  F.   Schaum's  Outline  of  Theory  and  Problems  of  Matrices. 
Schaum  Publishing  Co.,  New  York,  October  1962. 

12.  Lehmer,  D.  N.   A  Machine  Method  for  Solving  Polynomial  Equations. 
Journal  of  the  Association  for  Computing  Machinery,  Vol.  8,  No.  2, 
April  1961. 

13.  NASA  Progress  Report,  1965.   Prepared  under  Contract  No.  2-2609  by 
the  University  of  Santa  Clara,  Santa  Clara,  California. 

14.  Kokotovic,  P.  and  Siljak,  D.  D.   The  Sensitivity  Problem  in 
Continuous  and  Sampled  Data  Linear  Control  Systems  by  Generalized 
Mitrovic  Method.   IEEE  Transactions  on  Applications  and  Industry, 
Vol.  83,  No.  74,  September  1964,  pp.  321-324. 


71 


APPENDIX  I 


TABLE  OF  THE  CHEBYSHEV  FUNCTIONS  OF  THE  SECOND  KIND 


The  Chebyshev  functions  of  the  second  kind,  U ,  (C) »  are  defined 


by  the  recursion  relation: 


W5)  " 2  V?)  +  \-i(?>  ■  ° 


where: 


uQ(c)  =  o 
u1(c)  =  1 


These  functions  for  k  =  -1 ,  10  are 


k 

-1 
0 
1 
2 
3 
4 
5 
6 
7 
8 
9 

10 


ukU) 


-1 

0 

1 
2 


4c"  -  1 

8?3  -  4c 

4    2 
16c   12?  +  1 

32c5  -  32c3  +  6c 


64c 
128c' 
256c( 


512c' 


-  80c4  +  24c2  -  1 

-  192c5  +  80c3  -  8c 

-  448c6  +  240c4  -  40c2  +  1 

-  1024c7  +  672c5  -  160c3  +  10c 
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APPENDIX  III 


EXPANSION  OF 


B1      B2 
°1     °2 


Given  a  characteristic  polynomial  of  the  form: 


f(s)  *    y~    (bko<+  c^  +  d^)  s1 

k*o 


(A-1) 


Develop  a  general  expression  for: 


B2     C2 


'    B1C2  "  B2°1 


(A-2) 


where : 


n  n 

b,  -21  (-1  )\  ^  _,  B2 '  2_  (-1  >\  «> 

k»o  k»o 


00 


(a-3) 


k*o 


n 


c2.^(-,)\a>(k) 

k*o 


Consider  n       2: 


Equation  (A-1)  becomes: 


(b20(+c2p*i2)s2  +  (bl  01  +0^+^)8  +  (b^cft^)     «    o 
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Evaluating  equations  (A-3)  yields  : 
2 


B,  =  -to     +  b0  cj 

o  2     n 


B, 


«  -b,  cj     f  b„cj  2(2  X) 
1      n         2     n 


Cr-Co  +  C2"n 


L       -c,  o      +  c„  cj     (2  5  ) 

2  1      n         2     n 


From  which: 


B1C2-B2C1     "    Vlfln-VlUnSC2'l»?(2r)-b1C.fl^ 

b2coUn2(2l!)  +  b1  °2  "n3     *      ° 


After  grouping  terms,    this  is    : 


V2-B2C,    " 


'Vrbi0o)4,a->boWo)"n(2SJ 


(A-4) 


Consider  n     3   : 


Equation  (A-1 )  gives  : 


(b5*+c3p+d5)s3  +  (o20<+c2p +d2)s2  +  (^  oC+c^  td^s  +  (bo^+cop+do)»  0 


Evaluating  equations  (^-3)  yields  : 

B4  ■  -b     f  b0co  2  -  b,  co  5(2^S) 
1  o  2     n  3     n  ' 


B2  *  "b1  °n  *  b2  "n^2*  )  "  b3  "n^4  S  ^ } 


C     ■   -c     +  c_  &  2  -  c,  o  5(2*) 
o         2     n  3     n 


V    -c1u„  +  c^„2<2j)-c3U/(n2-,) 
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i\iote  that  in  each  case  only  one  term  of  the  summation  need  be  calculated, 
i.e; 


1(n)      V-1J         n  n  (n"° 


(A-5) 


Multiplying  as  indicated  in  equation  (A-3)»  cancelling  and  grouping 
terms  gives  : 


B,C2-B20, 


3/ ...2 


♦(»1o2-*201)*,n'"  (b1°3"Vl)<Un4<2,5)         ° 
♦  (b2c?-b5c2)  «n5         0  0 


Note  that  only  the  terms  uivolving  b-  and  c,  are  new,  that  is  : 

<Bi°2-B2°i>(n)-  ("i0a-Vi)(|lpl)*^5W.>  ^i^A1  ■T>>*) 


Consider  n  4  : 


evaluating  equations  (A-1 )  in  accordance  with  (A-5)  yields  : 

B  »     B     f  b  «,  4(4**-1) 

1  ^nO)   4  n 

B  .  B     ♦  b  «| 4(8X3-4t) 

2  2(nO)   4  " 

1        V?)    4  n 

C  *   C     *c  «4(8X3-4|) 
2     2(n-3)  4  n 
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Substituting  in  equation  (A-2),  cancelling  and  grouping  terms  gives  : 


*.w,  ■ 


Kxtension  to  the  general  case  : 


[kvM.]".. 


♦  [^-k|ej(-.fuj(-i)^  ♦ +&vfc>T«lfrfc" 


».VM,  • 


♦IvAvk"" • *kv  iJ,v.]Mi"«Jw^'-1 


L 


+  (\,c  -t  c  "Id 


a*- 1 


where   i   1,2, n   is  the  row  number 


j   1,2, ,n   is  the  column  number 


U-6) 
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Or,  in  compact  notation   : 


B,02-32C, 


ce 


i.J 


id 


U-7) 


where : 


CE    * 


ce. .     is  te&  (n  x  m)  upper  triangular  matrix 
whose  elements  are    : 


fi-ici+j-i  -bi+j-ici-i]  (-d^tuV)    n*+** 


i+j  <  a-2 
itj  £  n-2 


U-6; 
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APPENDIX  IV 


PROGRAM  PARAM  S 
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COMPUTED  DATA  FOR  SINGULAR  LINES 
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